Charge Disproportionation and Spin Ordering Tendencies in Na^Co-O^ 



O 

o 

(N 



m 

(N 



c3 



K.-W. Lee, 1 J. Kunes, 1 ' 2 and W. E. Pickett 1 

J Department of Physics, University of California, Davis CA 95616 
Institute of Physics, Academy of Sciences, Cukrovarnickd 10, CZ-162 53 Prague, Czech Republic 

(Dated: February 2, 2008) 

The strength and effect of Coulomb correlations in the (superconducting when hydrated) x ~l/3 
and "enhanced" x «2/3 regimes of Na^CoCb are evaluated using the correlated band theory LDA+U 
(local density approximation of Hubbard U) method. Our results, neglecting quantum fluctuations, 
are: (1) allowing only ferromagnetic order, there is a critical U c = 3 eV, above which charge dispro- 
portionation occurs for both :r=l/3 and £=2/3, (2) allowing antiferromagnetic order at £=1/3, U c 
drops to 1 eV for disproportionation, (3) disproportionation and gap opening occur simultaneously, 
(4) in a Co 3+ -Co 4+ ordered state, antiferromagnetic coupling is favored over ferromagnetic, while 
below Uc ferromagnetism is favored. Comparison of the calculated Fermi level density of states 
compared to reported linear specific heat coefficients indicates enhancement of the order of five for 
x ~0.7, but negligible enhancement for x ~0.3. This trend is consistent with strong magnetic be- 
havior and local moments (Curie- Weiss susceptibility) for x >0.5 while there no magnetic behavior 
or local moments reported for x <0.5. We suggest that the phase diagram is characterized by a 
crossover from effective single-band character with U » W for x > 0.5 into a three-band regime 
for x < 0.5, where U — > U e fj < U/v3 ~ W and correlation effects are substantially reduced. 

PACS numbers: 71.28.+d,71.27.+a,75.25.+z 
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I. BACKGROUND 

Since the discovery of high temperature supercon- 
ductivity in cuprates, there has been intense interest 
in transition metal oxides with strongly layered (quasi) 
two-dimensional (2D) crystal structures and electronic 
properties. For a few years now alkali-metal interca- 
lated layered cobaltates, particularly Naa;Co02 (NxCO) 
with x ~ 0.50 — 0.75, have been pursued for their 
thermoelectric propertiesJQ_The recent discovery^ and 
confirmation^, H % H 0, G3 EH of superconductiv- 
ity in the system NazCoOVyH^O for x sa 0.3 when inter- 
calated with water at the y w 0.3 level, has hei ghte ned 
interest in the NxCO system. The structure [T2T 1 1 3t IT^j 
is based on a 2D C0O2 layer in which edge-sharing C0O6 
octahedra lead to a triangular lattice of Co ions. Na do- 
nates its electron to the C0O2 layer, hence x controls 
the doping level of the layer: x=0 corresponds to Co 4+ , 
S=i low spin ions with one minority tig hole, and x — 1 
corresponds to non-magnetic Co 3+ . 

Nearly all reports of non-stoichiometric materials 
quote values of x in the 0.3 - 0.8 range, and the ma- 
terials seem generally to show metallic conductivity. 
The x=l endpoint member NaCo02, with rhombohe- 
dral R3m spacegroup lTol llr l is reported to be a con- 
ventional semiconductor . [1 71 1 1 8| The isovalent compound 
LiCo02 has been more thoroughly studied, with the con- 
clusion that it is a nonmagnetic band insulator with 2.7 
eV bandgap.fl^ l2fj| The x — endpoint has been an- 
ticipated by many to be a t\ g Mott insulator but is less 
studied; in fact, the Co 4+ formal oxidation state in a 
stoichiometric compound is rare. The sulfur counterpart 
C0S2 is metallic itinerant ferromagnet, close to being half 
metallic. A decade ago, Tarascon and coworkers reported 
synthesis of C0O2 but were unable to identify a complete 



structure. They concluded initially that the Co ions lay 
on a distorted triangular lattice [2U [2^]. More recently, 
further study by Tarascon et a/.|23j has traced the diffi- 
culty in pinning down the structure to the existence of 
two phases of C0O2, one stoichiometric and the other 
having 4% oxygen vacancies. C0O2 samples are metallic 
and nonmagnetic, hence cannot be said to contain Co + 
ions.01 

Much has been made of the similarities and differences 
of this cobaltate system compared to the cuprates. Both 
are layered transition metal oxide materials whose con- 
ductivity is strongly anisotropic. Both are in the vicinity 
of a Mott insulator (although the cobaltate one - C0O2 - 
is not well characterized). It is possible to vary the carrier 
concentration (x in the cobaltate formula) in both sys- 
tems, with the range in the cobaltates yet to be agreed 
on. In both systems there are specific superconducting 
regions: in the cuprates it is a "dome" 0.06 < x < 0.22, 
roughly, while in the cobaltates there are reports both of 
a dome 0.27 < x < 0.33 [13 and of a T c =4.5 K plateau 
for 0.28 < x < 0.37.[2^| However, the differences between 
the cobaltates and cuprates are substantial and expected 
to be crucial. Cobalt forms a triangular lattice, which 
frustrates antiferromagnetic (AFM) order, while the bi- 
partite square Cu lattice invites AFM order. The Co0 6 
octahedra are edge-sharing, rather than corner-sharing, 
making the bandwidth much narrower and the exchange 
coupling smaller. The cobaltates are electron-doped from 
the (anticipated) Mott insulator, as opposed to the most 
common hole-doped cuprates. And most striking, possi- 
bly: in the cobaltates T™ a:c =4.5 K compared to T™ ax = 
130 K (or higher under pressure) in cuprates. 

Another system for comparison is the transition metal 
disulfide based one, with Na 1/ / 3 TaS2-;t/H20 as the pri- 
mary comparison. In the (Nb,Ta)(S,Se)2 system, charge- 



density waves compete with superconducting pairing for 
the Fermi surface, with coexistence occurring in certain 
cases. The structure of the (for example) TaS2 layer is 
identical to that of the C0O2 layer, consisting of edge- 
sharing transition metal chalcogenide octahedrons. In 
these dichalocogenides, as in the cobaltate system, two 
well defined stages of hydration have been identified. j2|| 
In the first stage H2O is incorporated into the same layer 
as the cation (typically an alkali ion), and in the second 
stage two H2O layers are formed on either side of the 
cation layer. The similarity in increase in the c lattice 
parameter compared to the cobaltates is illustrated in 
Fig. ID 

The electron concentration in Na 1 / 3 TaS2-yH20 is spec- 
ified by the Na concentration £=|, and in this system 
superconductivity occurs at 4-5 K (the same range as 
in the cobaltates) regardless of the concentration y of 
water molecules intercalated into the structure. [27] Spe- 
cific thermodynamically stable phases were identified at 
y = 0, §, 0.8, 1.5, and 2.(23, E[J The level of elec- 
tron donation seems to be crucial: using Yjyg and 
based on the trivalcnt ions leads to the same supercon- 
ducting transition temperature. Using the divalent ion 
Mn (d 5 ,S = §), at the same electron donation level 
Mn 1 / 6 TaS2 is ferromagnetic (FM). Intercalating this FM 
compound with H2O leads again to T c ~ 4 K. This latter 
behavior is understood as the water-induced separation 
of TaS2 layers decreasing the interlayer magnetic cou- 
pling sufficiently to inhibit long-range magnetic order, 
thereby allowing the innate superconducting tendency in 
the doped TaS2 layer to assert itself. 

Much of the emphasis, both experimental and theoret- 
ical, has been directed toward the superconducting be- 
havior of the cobaltates, but the long- known behavior of 
the tantalum disulfides just mentioned suggests the su- 
perconductivity may not be so distinctive. Reports of 
the magnetic behavior in these cobaltates have been of 
particular interest to us. Except for a charge dispro- 
portionated and charge-ordered phase in a narrow range 
around x=0.5[3l| identified by its insulating behavior, all 
samples are reported to be metallic. For x in the 0.5-0.8 
range, the susceptibility x(T) is Curie- Weiss-like with re- 
ported moment of order 1 /j,b per Co 4+ . 0. [J- E} This 
local moment is normally interpreted as indicating the 
presence of correlated electron behavior on the Co sub- 
lattice, and most theoretical treatments have assumed 
this viewpoint. 

Some phase transitions have been reported in the high 
x region. Magnetic ordering at 22 K with small or- 
dered moment has been reported for ir=0.75j3Jj based 
on transport and thermodynamic data, and the same 
conclusion was reached by Sugiyama et al. from /xSR 
studies. 0j H3 Boothroyd et al. performed inelastic neu- 
tron scattering on x=0.75 single crystals and observed 
FM spin fluctuations. [3(| Field dependence of the ther- 
mopower measured by Wang et al. 10] indicated that the 
spin entropy of the magnetic Co system {i.e. the spins 
of the Co 4+ ions) is responsible for the unusual thermo- 
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FIG. 1: Change in the c axis lattice parameter with ad- 
dition of H2O, in the two systems Nao.aCoCvH^O and 
Na 1/ / 3 TaS2-H20, illustrating the great similarity. For cobal- 
tates, data are from Foo et al. l32j ( empty circles), Jin et aZ.0 
(filled circle), and Schaak et aZ.[Ol (asterisk). For the chalco- 
genides, date are from Johnston J23] (empty diamonds) and et 
al.fx% ( filled diamonds'). 



electric behavior. For a;=0.55, Ando et al. reported 37] a 
rather large linear specific heat coefficient 7=56 mJ/mol- 
K 3 . Thus for x > 0.5 magnetic Co ions are evident and 
are strongly influencing the electronic low energy excita- 
tions. 

However, for samples with x ~0.3 (i.e. the supercon- 
ducting phase) many reports concur that the Curie- Weiss 
behavior of \ vanishes. 0, H, El US Ell] • In addition, the 
specific heat 7 is much smaller, with values around 12 
mJ/mol-K 2 reported. 0, 0] It is extremely curious that 
local moments should vanish so near to what has been be- 
lieved would correspond to a Mott insulator (x=0, Co 4+ 
in C0O2), and that superconductivity appears only in the 
moment-free regime. In the strongly interacting single- 
band triangular lattice picture, the x = system corre- 
sponds to the half-filled triangular lattice that has been 
studied extensively for local singlet (resonating valence 
bond) behavior. |40j The ground state of that model has 
however been found to be Neel ordered. ^lj In any case, 
the x « 0.3 region of superconductivity in NxCO is how- 
ever well away from the expected Mott-insulating regime, 
and the behavior in such systems is expected to vary 
strongly with doping level. 

Much of the language used in characterizing this sys- 
tem (above, and elsewhere) has been based on the local 
orbital, single band picture. As discussed more fully be- 
low, the doping in this system occurs within the threefold 
tig complex of the Co ion, with degeneracy only slightly 
lifted by the non-cubic structure. The question of single- 
band versus multiband nature of this cobaltate system 
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is possibly one of the more important issues to address, 
because it can affect strongly the tendency toward corre- 
lated electron behavior. 

Although the primary interest has been in the super- 
conductivity of NxCO, there is first a real need to under- 
stand the electronic structure of the normal state of the 
unhydrated material, and its dependence on the doping 
level x. The electronic structure of the x=l/2 (with or- 
dered Na) compound in the local density approximation 
(LDA) has been described by Singh. [H|43| Within LDA 
all Co ions are identical ("Co 3 5+ "), the Co ti g states are 
crystal-field split (by 2.5 eV) from the e g states, and the 
tig bands are partially filled, consequently the system is 
metallic consistent with the observed conductivity. The 
tig band complex is W ~ 1.5 eV wide, and is separated 
from the 5 eV wide O 2p bands that lie just below the 
Co d bands. Singh suggested that the expected on-site 
Coulomb repulsion (which has not been calculated) is 
[7=5-8 eV on Co, which gives U » W so that correla- 
tion effects can be anticipated. 

Notwithstanding the experimental evidence for non- 
magnetic Co ions in the superconducting material, most 
of the theoretical discussion H3, |H |iE EE |H |4^, E3, E3 
has focused on the strongly interacting limit, where U 
is not only important, but in fact is presumed to be so 
large that it prohibits double occupancy, as described by 
the single band t — J model. The lack of local moments 
and only weak to moderate enhancement of the specific 
heat 7 suggests that a more realistic picture may be re- 
quired. Undoubtedly the single band scenario is a limited 
one: although the rhombohedral symmetry of the Co site 
splits the t 2g states into a g and e' g representations, the 
near-octahedral symmetry makes them quasi-degenerate. 
Koshibae and Maekawa have shown that the band dis- 
persion in the ti g band complex in these cobaltates dis- 
plays unexpected intricacies, including some analogies to 
a Kagome lattice. £3 

In this paper we begin to address the correlation 
question by taking the strongly correlated viewpoint 
and using the correlated band theory LDA+Hubbard 
U (LDA+U) method. We investigate two distinct re- 
gions of the phase diagram by focusing on x =1/3, the 
regime where superconductivity emerges, and x =2/3 
where more magnetic behavior is observed. We find that 
U > U c = 3 eV leads to charge disproportionation (CD) 
and gap opening for both £=1/3 and x=2/3 if only FM 
order is allowed. For the Neel ordered case at £=1/3, 
the corresponding transition occurs at U c = 1 eV. The 
availability of three distinct sublattices for the ordering, 
coupled with strong 2D fluctuations, may destroy long- 
range order and make local probes important in studying 
charge disproportionation and correlation. The small val- 
ues of U c that we obtain even for x—1/3 tend to confuse 
the theoretical picture, since there seems to be a con- 
spicuous absence of correlated electron behavior in this 
regime of doping. 




percell. The unconnected large spheres represent nonmag- 
netic Co 3+ ions, while the large and small connected spheres 
represent oppositely directed Co 4+ spins when ordered anti- 
ferromagnetically. 



II. STRUCTURE AND METHOD OF 
CALCULATION 

Our calculations are based on the hexagonal structure 
(space group P6322), obtained by Jansen and Hoppe, 
having lattice constants (a = 2.84A, c = 10.81 AL|53j The 
supercell (\/3a x \/3a x c/2) is used so that at the concen- 
tration x — 1 that we consider, two (or possibly three) 
inequivalent Co ions, viz. Co 3+ and Co 4+ , are allowed 
to emerge in the process of self-consistency. The allowed 
order is displayed in Fig. [21 Since we are not analyz- 
ing the very small intcrlayer coupling here, a single layer 
cell is used. In the supercell (space group PSlm, No. 
157), atomic coordinations are Na at the la (0, 0, 1/2) 
above/below the Co site at the la (0, 0, 0), and the other 
Co sites are the 2b (1/3, 2/3, 0). Oxygen sites are the 3c 
(2/3, 0, z ) and the 3c (1/3, 0, z ) positions, respectively. 
The O height z = 0.168(c/2) = 0.908A, which is relaxed 
by LDA calculation, [42] produces the Co-O-Co bond an- 
gle 98.5° (90° for undistorted) , so that the octahedra is 
considerably distorted. 

Two all-electron full-potential electronic methods have 
been used. The full-potential linearized augmentcd- 
plane-waves (FLAPW) as implemented in Wien2k code 
|54j and its LDA+U [55l l5t| extension were used. The 
s,p, and d states were treated using the augmented 
plane wave+local orbitals (APW+lo) scheme [57j , while 
the standard LAPW expansion was used for higher Vs. 
Local orbitals were added to describe Co 3c? and O 
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FIG. 4: Band structure (in the virtual crystal approximation) 
along high symmetry lines (left panel) and the aligned density 
of states (right panel) for the x — i cobaltate in the local 
density approximation. The a g symmetry band is emphasized 
with circles proportional to the amount of a g character. The 
a g density of states is indicated by the darker line. 



FIG. 3: Local density approximation bands of Na^OoO^ in 
the virtual crystal approximation where there is one Co ion 
per cell, shown along the principle high symmetry directions. 
The e g bands lie above 1.5 eV; the bands below -1.5 eV are 
predominantly oxygen 2p in character. The thickened lines 
emphasize the bands within the ti g complex (-1.5 to 0.5 eV) 
with strong a g character. 



2s and 2p states. The basis size was determined by 
RmtKmax = 7.0. The full-potential nonorthogonal local- 
orbital minimum-basis scheme (FPLO)[58|, EH was also 
used. Valence orbitals included Na 2s2p3s3p3d, Co 
3s3p4s4p3d, and O 2s2p3s3p3d. The Brillouin zone was 
sampled with regular mesh containing 50 irreducible k- 
points. Both popular forms [60l l6ll| of the LDA+U func- 
tional have been used in our calculations, with no impor- 
tant differences being noticed. 



III. WEAKLY CORRELATED LIMIT. 

LDA electronic structure at x = | . The e g — t 2g crys- 
tal field splitting of 2.5 eV that can be seen in the full 
band plot in Fig. places the (unoccupied) e g states 
(1 eV wide) well out of consideration for most low en- 
ergy effects. The O 2p band complex begins just be- 
low the bottom of the t 2g bands (see Fig. and is 5.5 
eV wide. The states into which holes are doped from 
NaCo02 come from the 1.6 eV wide t2 g band complex. 
The trigonal symmetry of the Co site in this triangular 
C0O2 layer splits the t2 g states into one of a g symmetry 
[{\xy > + \yz > + \zx >)/v / 3 in the local CoOg octahe- 
dron coordinate system] and a doubly degenerate pair e' g 
[{\xy > +a\yz > +a 2 \zx >)V3 and its complex conju- 



gate, where a — exp(2iri/3)]. 

The t 2g band complex that is intersected by the Fermi 
level is shown in more detail in Fig. ^ where the 
bands with primarily a g character are shown in the "fat- 
bands" representation against the corresponding densi- 
ties of states. The band dispersions agree well with those 
calculated by Rosner et aL[63| The a g character is strong 
at the bottom of the t 2g complex as well as at the top, 
and illustrates that holes doped into the band-insulating 
NaCo02 (x—1) phase go initially into one a g band that 
is rather flat for ^25-30% of the distance to the zone 
boundary. Based on a rigid band interpretation using 
this £=1/3 density of states (DOS), doped holes enter 
only a g states until x w0.6, whereupon an e' g Fermi sur- 
face begins to form. This observation is consistent with 
the ir=0.5 Fermi surface shown by Singh[42l] which has 
six e' g cylinders. 

It is of interest to view this band structure from the 
viewpoint of a single isolated tight-binding s-band on a 
triangular lattice with near neighbor hopping only, which 
is intended[il El El |H |H lEH to model the a g band 
dispersion. The a g DOS lies higher than that of e' g not 
because its band center lies higher (in fact its centroid 
is somewhat lower) but rather due to the particular dis- 
persion and to a substantially larger effective bandwidth. 
Judged from the dispersion curves themselves, the a g and 
e' g bands differ little in width. However, nearly all of the 
e' g states lie within a 1.0 eV region, whereas the a g DOS 
extends over 1.5 eV. 

The a g band dispersion in Fig. 0] does resemble that 
of the simple tight-binding model — t X]<ij> ( c l c j + h.c.) 
with a negative value of t. (The band structure also indi- 
cates that the e' g hopping integral has the opposite sign 
to that of a g .) The a g -projected DOS however is noth- 
ing like that of the tight-binding model. 0] The reason 
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is twofold. First, there is mixing of the a g with the e' g 
bands over most of the Brillouin zone. The hybridization 
is evident along the M-r line in Fig. it is less obvi- 
ous along the T-K line because the mixing happens to be 
accidentally small for the 2;= 1/3 C0O2 layer structure. 
For other Co-0 distances and bond angles, and for £=0 
doping level (not shown) , the mixing of the a g band with 
the lower e' g band becomes much larger. A second reason 
for the actual shape of the DOS is due to the influence 
of some second- neighbor hopping, which makes the 
a g band near fc=0 much flatter than the tight-binding 
model, or even disperse slightly upward before turning 
downward. 

Some details of the band structure should be clarified. 
The upward dispersion of the a g band around the T point 
(mentioned above) also seems to be affected by interlayer 
coupling, which can depress the band at k=0. Johannes 
and Singh have reported that, even for well separated 
C0O2 layers (i.e. when hydrated) the a g band may still 
disperse upward 63j before turning down. Even for C0O2 
layer geometries for which there is no upward curvature, 
the a g band remains unusually flat out to almost 1/3 of 
the way to the zone boundary. Either behavior is indica- 
tive of extended hopping processes. 

Magnetic Order with LDA. Analogous to the results of 
Singh for x — 0.3,0.5,0.7 0jH3|, we find ferromagnetic 
(FM) tendencies for £=1/3 within LDA. In disagreement 
with experiment (no magnetic order is observed for x ~ 
0.3, nor even any local moment at high temperature) 
a half metallic FM result is found, with a moment of 
2 hb I supercell that is distributed almost evenly on the 
three Co ions. The exchange splitting of the ti g states 
is 0.5 eV, and the Fermi level (Ep) lies just above the 
top of the fully occupied majority bands (the minority 
bands are metallic). The FM energy gain is about 45 
meV/Co. With the majority bands filled, the filling of 
the minority ti g bands becomes |, leading to larger e' g 
hole occupation than for the paramagnetic phase. Hence, 
unlike the standard assumption being made so far, x = i 
is a multiband (a g + e' ) system (within LDA, whether 
ferromagnetic or paramagnetic). Attempts using LDA 
to obtain self-consistent charge disproportionated states, 
or AFM spin ordering, always converged to the FM or 
nonmagnetic solution. 

Fermiology. Suspecting from the 5=1/2 spins and the 
two-dimensionality that quantum fluctuation is an im- 
portant aspect of this system, it is possible that the 
x—1/3 system is a "fluctuation-induced paramagnet" due 
to the lack of account of fluctuations in the electronic 
structure calculations. Whatever the underlying reason, 
a page can be taken from the high T c cuprate chapter 
of materials physics that, even in the presence of consid- 
erable correlation effects, in the magnetically-disordered 
metallic phase the paramagnetic Fermi surface (FS) will 
emerge. The lack of any observed magnetic behavior in 
the x ~0.3 region reinforces this expectation. 

In Fig. show the x—1/3 LDA FS, which is sim- 

ilar to the x = 0.5 one shown by Singh 42]. A large 




FIG. 5: Fermi surface for (virtual crystal) NxCO, a;=0.30, in 
the two dimensional Brillouin zone. The large cylinder con- 
tains a g holes, whereas the six small cylinders contain holes 
that are primarily eg-like. 

T-centered hole cylinder (mean radius Kp) shows some 
flattening perpendicular to the T-K direction, this cylin- 
der holds 0.43 a g holes/Co. In addition, there are six 
additional, primarily e' g in character, hole cylinders lying 
along the T-K directions, containing 0.04 holes in each 
of the six small cylinders (radius kp). The total is the 
0.67 holes necessary to account for the x = 0.33 elec- 
tron count. This FS geometry leads to several important 
phase space features. There are the nesting wavevectors 
that translate one of the small cylinders into another, 
giving three distinct intercylinder nesting vectors as il- 
lustrated in Fig. [S] If these cylinders were circular, these 
vectors would represent strong nesting vectors for charge- 
or spin-density waves. In addition, the susceptibility for 
Q < 2kp intra surface scattering processes is constant in 
two dimensions. [67j The calculated cylinders have an ec- 
centricity of 1.25, weakening these nesting features some- 
what. There are in addition the corresponding processes 
with Q < 2Kp of the large cylinder. 



IV. INCLUSION OF CORRELATION EFFECTS 

Despite the feature of the LDA+U method that it 
drives local orbital occupations to integral occupancy (as 
U increases) , to our knowledge it has never been used to 
study the question of charge disproportionation. In this 
section we show that moderate values of U lead to CD 
at both a;=l/3 and x—2/3. For the triangular lattice, 
threefold expanded supercells (-s/3 x -s/3, see Sec. II) are 
convenient, and 3;= 1/3 lies very close to the supercon- 
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FIG. 6: Effect of the intraatomic repulsion U on the magnetic 
moment (left axis) and energy gap (right axis) of the Col and 
Co2 ions in the supercell. Top panel: result for ferromagnetic 
order. Changes in the magnetic moments, indicating dispro- 
portionation to formal charge states Co 3+ and Co 4+ , begins 
at (7 C = 3 eV. The opening of the gap (half metallic to in- 
sulating) occurs simultaneously at U c . Bottom panel: results 
for antiferromagnetic order. Already for U c =l eV, the Co 4+ 
moment is large (the Co 3+ moments is zero by symmetry) 
whereas the gap is just beginning to open. 



ducting composition while x=2/3 is representative of the 
x > 0.5 region that shows correlated behavior. [6^ 

LDA + U magnetic structure and energies. The behav- 
ior of the LDA+U results (Co moments and the energy 
gap) versus U was first studied for x=l/3 (on-site ex- 
change was kept fixed at the conventional value of 1 eV 
as U was varied). The dependence of the magnetic mo- 
ment and band gap on U for FM ordering is shown in the 
top panel of Fig. For U < U c = 3 eV, the moments on 
the two inequivalent Co sites are nearly equal and sim- 
ilar to LDA values (which is the U — >0 limit). Above 
U c , disproportionation into S = h Co 4+ and 5 = Co 3+ 




U(eV) 

FIG. 7: Effect of the intraatomic repulsion U on the mag- 
netic moment (left axis) and energy gap (right axis) of the 
Col and Co2 ions in the supercell for x=2/3. Changes in the 
magnetic moments, indicating disproportionation to formal 
charge states Co 3+ and Co 4+ , begins at U c = 3 eV just as 
in the corresponding a;=l/3 case. Note that applying small 
but increasing U decreases the moment somewhat until dis- 
proportionation occurs. The opening of the gap (half metallic 
to insulating) occurs simultaneously at U c . 



ions is nearly complete at U = 3.5 eV and is accompanied 
by a metal-insulator (Mott-like) transition from conduct- 
ing to insulating. The gap increases linearly at the rate 
dE g /dU = 0.6 for U > 3.5 eV. For the (insulating) [7=5 
eV case, nonmagnetic Co 3+ states lie at the bottom of 
the 1.3 eV wide gap, with the occupied Co 4+ e' g states 1-2 
eV lower. The spin-half "hole" on the Co 4+ ion occupies 
the a g orbital as expected. 

In our choice of (small) supercell, this CD is necessar- 
ily accompanied by charge order, resulting in a honey- 
comb lattice of spin half ions. In a crystal there would 
be three distinct choices of the ordered sublattice (corre- 
sponding to the three possible sites for Co 3+ ). Of course, 
even at a rational concentration such as £=1/3, CD may 
occur without necessarily being accompanied by charge 
ordering when thermal and quantum fluctuations are ac- 
counted for. Regarding the disproportionation, we note 
that, based on the Mullikan charge decomposition in the 
FPLO method, the charges on the "Co 3+ " and "Co 4+ " 
ions differ by only 0.25-0.3 electrons. This small value 
reflects the well known result that the formal charge des- 
ignation, while being very informative of the magnetic 
state and indispensable for physical understanding, does 
not represent actual charge accurately. 

The analogous calculation can be carried out allowing 
AFM order of the Co 4+ ions, and the results are shown 
in the bottom panel of Fig. A Co 4+ moment grows 
(disproportionation) immediately as U is increased from 
zero. The gap opens around U = 1 eV and increases at 
the rate dE g /dU = 0.4. Thus for AFM spin order, the 
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FIG. 8: Evolution of the Co 3+ and Co 4+ 3d states as charge 
disproportionation occurs, for x=l/2> and FM spin order (top 
panel) and AFM spin order (bottom panel). Note that the 
evolution in the FM case is from a half metallic FM. 



critical value is no higher than U c = 1 eV. 

At 2=2/3, CD will lead to only one Co 4+ ion in the 
supercell, so only FM ordering can be considered. The 
corresponding behavior of the moment and the gap are 
presented in Fig. The Co moments remain nearly 
equal but slowly decrease from their LDA value up to 
U c = 3 eV, whereupon again CD occurs abruptly. The 
moments are "well formed" by U — 4 eV but continue 
to evolve somewhat beyond that. In this case dE g /dU = 
0.67. 

Our LDA+U results, showing charge disproportiona- 
tion for U c = 3 eV (FM) or U c = 1 eV (AFM) are very 



different from earlier reports where little change was ob- 
tained even for larger values of U when symmetry break- 
ing by disproportionation was not allowed. |64j This dif- 
ference serves as an alert that LDA+U results can be 
sensitive to what degrees of freedom are allowed. 

It is far from obvious that charge disproportiona- 
tion and gap opening should occur simultaneously with 
LDA+U, although physically the phenomena are ex- 
pected to be closely related. The evolution of the Co 3+ 
and Co 4+ 3d states with U in the critical region is shown 
in Fig. [S] In the FM case, the system evolves from a 
half metallic configuration, still visible for U = 3 eV, and 
gap opening occurs just when the minority a g band con- 
taining 2/3 hole per Co ion splits off from the valence 
band. At the point of separation, the a g states split into 
an unoccupied narrow band containing one hole for each 
of the two Co 4+ ions, and an occupied (and also narrow) 
band on the Co 3+ ion. This disproportionation can be 
identified from the strong change in the occupied states 
around -0.6 to -0.3 eV. Thus while disproportionation in 
principle could occur before the gap opens, it does not 
do so here. 

When symmetry is broken according to the AFM or- 
dered (and disproportionated) superlattice shown in Fig. 
13 the critical value of U is 1 eV or less. This "ease" in 
gap opening is no doubt encouraged by the narrow band- 
width of the unoccupied a g band. A spin up Co 4+ ion 
is surrounded by three spin down Co 4+ ions and three 
Co 3+ ions, neither of which have a g states of the same 
spin direction at the same energy. The surviving band- 
width reflects the effective coupling between a g states on 
second-neighbor Co ions. For this AFM ordered phase, 
the DOS of Fig. [S] indicates also little or no dispropor- 
tionation before the gap begins to open. 

Exchange Coupling. The Coulomb repulsion parame- 
ter U has not been calculated for these cobaltates, but 
several estimates for other cobaltates put U at 5 eV or 
above. ^3 EH E3 On the other hand, Chainani et al. fit 
cluster calculation results to xray photoemission data on 
samples showing both charge states, and found that U 
in the 3-5.5 eV range work equally well.|s3 This range 
is rather lower than what has generally supposed, thus 
the appropriate value of U is quite uncertain. Here we 
concentrate on [7=5.4 eV results, but our calculated be- 
havior is not sensitive to variation of U in this range. 
In this CD regime (also charge ordered, due to the con- 
straint of the supercell), AFM ordering gives 1.2 mRy/Co 
lower energy than does FM order. In terms of nearest 
neighbor coupling on the resulting honeycomb lattice, 
the FM - AFM energy difference corresponds to J = 11 
meV. Referring to the paramagnetic bandwidth identified 
above, the corresponding superexchange constant would 
be 4t 2 /U ~ 20 meV. Again, we note that the a g DOS dif- 
fers greatly from the single band picture that was used 
to obtain this value of i=0.16 eV. 
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V. DISCUSSION OF INTERACTION 
STRENGTH 

Our calculations indicate that, as long as U > 3 eV as 
is generally thought, at both x=l/3 and x=2/3 there is a 
strong tendency to disproportionate, with one result be- 
ing a Co 4+ ion with a local moment. Disproportionation 
into an AFM honeycomb lattice occurs already by U c = 1 
eV in our mean field treatment. At least in the presence 
of charge order in a honeycomb arrangement, there is an 
AFM nearest neighbor exchange coupling J « 11 meV. 
The Neel state is known to be the ground state of the 
AFM Heisenberg honeycomb lattice. These charge order- 
ing tendencies may be expected to be strongly opposed 
by thermal and quantum fluctuations that are expected 
for low coordination and small spins in 2D layers. To 
date, disproportionation and charge ordering have only 
been reported for x=l /2.|3l) 

In addition to producing local magnetic moments, 
charge disproportionation might be expected to intro- 
duce new coupling to the lattice. Since the radii of Co 2+ 
and Co 3+ differ by 15% (0.74 A vs. 0.63 A), dispropor- 
tionation into those charge states would be expected to 
couple strongly to local oxygen modes. In octahedral co- 
ordination, however, the Co 4+ ion radius is almost indis- 
tinguishable from that of the Co 3+ ion,[6(| so there may 
be little evidence in the lattice behavior even if Co 3+ - 
Co 4+ charge disproportionation occurs. 

In spite of the prevalent theoretical presumption that 
correlation effects may be playing an essential role in the 
superconductivity of this cobaltate system, the data seem 
to be suggesting otherwise. In the x > 0.5 regime, indeed 
local moments are evident in thermodynamic and trans- 
port data, spin fluctuations have been observed by neu- 
tron scattering, and the linear specific heat coefficients 
are large, 7=48-56 mJ/mol-K 2 . Comparing this value to 
our calculated band value, j — 10±2 mJ/mol-K 2 , leads 
to a factor of five enhancement due to dynamic corre- 
lation effects. Magnetic ordering around x ~ 0.75 also 
attests to substantial correlation effects. Our finding of 
disproportionation for U > 3 eV (in mean field) is consis- 
tent with the experimental information and a correlated 
electron picture. 

In the superconducting regime x rs 0.3, the emerging 
picture is quite the opposite. The specific heat coeffi- 
cient is ordinary, with the reported values 0, Q cluster- 
ing around 7= 12-13 mJ/mol-K 2 indistinguishable the 
band value 70 13 mJ/mol-K 2 . In addition, there is no 
local moment (Curie- Weiss) contribution to the suscepti- 
bility, and other evidence of enhanced properties is lack- 
ing (magnetic field dependence of the resistivity is small, 
for example). In short, evidence of substantial correla- 
tion effects due to the anticipated strong on-site Coulomb 
repulsion U is difficult to find for x <0.5. Moreover, as 
discussed in the Introduction, the x=0 endpoint is not a 
Mott insulator, but rather a nonmagnetic metal. |24| 

It is essential to begin to reconcile the microscopic 
model with observations. There are several indications 



that the simple single band picture is oversimplified, one 
of the most prominent being that there is no evidence 
that the a g state is significantly different in energy from 
the e' g states, i.e. the t2 g degeneracy is still essentially in 
place. Due to the form of dispersion in the C0O2 layer, it 
is still the case that holes doped into the band insulator 
NaCo02 go into the a g band, making it viable to use a 
single band model in the small 1 — x regime with a rather 
robust value of U, with a value of U ~ 3 eV possibly be- 
ing sufficient to account for correlated electron behavior 
(W « 1.5 eV). 

The x < 0.5 regime seems to require reanalysis. It is 
quite plausible, based both on the LDA band structure 
and the observed properties, that for x < 0.5, the system 
crosses over into a three-band regime where the full ti g 
complex of states comes into play. The multiband nature 
tends to mitigate correlated behavior in at least two ways. 
Firstly, doped carriers that go into a multiband complex 
may simply find a smaller "phase space" for approaching 
or entering the Mott-Hubbard insulating phase, as the 
carriers have more degrees of freedom. Not completely 
separate, perhaps, is the extensive study of Gunnarsson, 
Koch, and Martin, (63, lfi!j|7fj] that strongly suggests that 
in a multiband system of N bands, the effective repulsion 
strength becomes U e ^ = U/VW. For these cobaltates 
with carriers in the t2 g bands, N=3, and U c < 3 eV would 
become — > 3 eV/y/3 ~ W, and correlation effects 
diminish considerably. Secondly, screening will increase 
as hole doping occurs from the band insulator x=l, re- 
ducing - perhaps strongly - the intra-atomic repulsion U 
to a value near W.. 



VI. SUMMARY 

In this paper we have begun an analysis, coupled with 
close attention to the observed behavior, of the strength 
of correlation effects in this cobaltate system that super- 
conducts when hydrated. In this initial work, we have 
used the mean field LDA+U method to evaluate the 
effects of Hubbard- like interactions in NxCO, and find 
charge disproportionation and a Mott insulating state for 
Coulomb repulsion U — 3 eV or less, for both x = 1/3 and 
x = 2/3, when fluctuations are neglected. Ferromagnetic 
tendencies for small U evolve to nearest neighbor antifcr- 
romagnetic coupling J w 11 meV for U ~ 5 eV, at least 
if charge disproportionation occurs. The only insulating 
phase reported so far has been at £=1/2, with strong 
evidence |3l| that it is due to charge disproportionation 
and charge order (and probably magnetic order). 

The x = | LDA FS has been described, following the 
presumption (based on the observation of at most mod- 
erately correlated behavior) that 2D fluctuations will re- 
store the paramagnetic metallic state. There are strong 
indications however that strong interactions, clearly evi- 
dent for x > 1/2, have become muted in the regime where 
superconductivity appears. On the one hand, the elec- 
tronic structure and FS indicate that multiband effects 
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must be considered in this regime, which in itself will 
decrease the effective repulsion U. Independently, U will 
be decreased by screening as the system becomes increas- 
ingly metallic. On the experimental side, the behavior of 
both the magnetic susceptibility and the linear specific 
heat coefficient point to a lack of "enhanced" behavior 
for x « 0.3. 
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